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/define KXN_HX 1 
/define MAX HI 2001 
/define hxJEnoex 1 

/define MAX_IDX 4000 

/* #def ine TISSUB */ /* TUsu^emiivalent .poxy plastic */ 

£L#Xe C4 ' /* C4 plastlo explosive •/ 

J?5sXiSe RDX */ /* XUJX sheet explosive */ 

'/! J«3iJe W */ /* ««ter Gel explosive •/ 

{! {deJS. S?H •/ /• 40% tfynamits stick */ 

/* new way of determining low */ 

/define si .0247 

/define *2 .01492 

/define s3 .265 

/define 14 112.6 

/define £5 25,198 

/define a6 .6218 

/define 17 .265 

/* def ine substance parameters */ 

/ifdef WO 
/define ol 9.732 
/define o2 6.108 
/define o5 1.218 
/define XO .547 
/define XL .961 
/endif 

/ifdef BDX 

/define ol 9.732 

/define o2 6.108 

/define o5 1.218 

/define X0 .65 
/define XL .86 
/endif 

/ifdef C4 
/define cl 9.732 
/define 02 6.108 
/define eS 1.218 
/define K0 .6522 
/define XL .87 
/endif 

/ifdef DYN 
/define cl 570.46 
/define o2 4.352 
/define c5 .304 
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/define KO .522 
/define KL .765 
/endlf 

/ifdef TISSUE 
/dafina cl 9798 
/define o2 3.8837 
/define 08 0.993 
/define KO .685 
/define KL .825 
/endif 

double bh (double km) ; 
double bh (double km) 
< 

return(cl*pov((kmfc5),c2)); 

double Kref (double Hi,double Km, double kO) ; 
double Kref (double Hi , double Km, double kO) 

return ( ( (Hi+bh (Km) ) *kO*KL) / ( (bh (Km) *KL) + (Hi*kO) ) ) ; 

double alpha (double km) ; 
double alpha (double km) 

i 

return ( (a 1+ (E2*km)-(a2*a3) ) / (km*km) ) ; 

double beta (double km); 
double beta (double km) 

ratum(C«4+((i6-km)*(z5/(z6-t7))))/km) ; 

double newlov (double h, double km) ; 
double nevlov (double h, double km) 

return (h* (1/ (km+ (alpha (km) * (h/ (h+beta (km) )))))); 

} 

double find_Ka (double hi /double Kair, double kref) ; 
double find^Km (double hi, double Kair, double kref) 
{ 

/* find the Km that approximates the desired Kref given high val,kO */ 
int x,bitval; 

double lsbval,approx_kref ; 
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Isbval - 0.8; 
bitval « o; 

for (x-0;x<8;x++) 
{ 

bitval- (bitval«i) 1 1; 
isbval - Isbval/ (double) 2.0; 

approxjcref - (Kref(hi,( (double) .l+( (double) bitval*lsbval)) # Kair))/ 

if (approxjcref < kref ) 
bitval-bitvalfi (Oxf e) ; 

) 

return (( (double) bitval*lsbvai) +. 1) ; 

double findKta Lov (double hi, double low); 
double findKtt_Lov (double hi, double low) 

(nt'J^biSSl? ***** afiproxl,aatea given high val,ko */ 

double labval , approxJLow; 

Isbval - 0.8; 
bitval - 0; 

for (x-0/x<8/x++) 

bitval- (bitval«i) 1 1; 
Isbval - Isbval/ (double) a,0; 

approx_lov - (Low (hi, ( (double) • 1+ ( (double) bitval*lebval) ) ) ) ; 

If (approx low < low) 

bltval=bitval4(0xfe) ; 

> 

return ( ( (double) bitval*lsbval) +. 1) ; 



/* create the histogram */ 

for (hint - MIN_HI; hint < MAX_HI; hint +- HIJTODEX) 

h - (double)hint; /* G « t hi double valxje.*/ 

/* Set up the header values and the KIdx */ 
Hdr£Hi VALUE] - hint; 
KIdx «"o; 
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/* Gat the hi and lo kref */ 
hi kref ■ Kref(h, 0.29, kO) ; 
loTaref = Kref(h, 0.8, kO); 
k-To_kref; 

lastl e -100.O; 
diffl - 1000.0; 
vhila (k<hijcref) 

km-find Kn(h,ko,k) ; 
kr»Xref?h,km,k0) ; 
l-Low(h,ka) ; 

if (((l-lastl)<diffl)M(km>.29)) 
diffl - 1 - lastl; 

lastl - l; 

if (b>800.0) 
{ 

k-k*l.04; 
> alas 

k-l.01*k; /* 1% bins */ 



/* do it again, bat use diffl to find values V 
k-lo kref; 
knaflnd Km(h;ko,k) ; 
l»Lov(h7ka) ; 
f indl-(int) l; 



/* adjust diffl to a power of 2 */ 
tdiffl-0; 

vhila ((1 « (tdiffl+1)) <- (int) diffl) 
tdiffl++; 

knofindKn Lov(h, (double) f indl) ; 
k-Kref (h,km,k0) ; 

/* Save the minimum low and the eoale factor */ 
Hdr[KXH L0] - findl; 
Rdr[L0_SCALS] - tdiffl; 

vhila (k < hi Jcref ) 

km°findXA Low (h, (double) findl) ; 
k-Kraf(h,fcn,kO); 

/* Save the necessary information into the values */ 
Kref Tab [KIdx] «= (float) k; 
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KXdx++; 

/* increment low */ 
findl (l « tdlffl); 

/* increment bin count */ 
bincnt+-l; 

> 

/* Nov we have the table, write out the header then the table */ 
HdrfKAX LO] - findl; 

bwritten - write (fhndl, (char *)Kdr, sizeof (int) *4) ; 
if (bwritten I- (sizeof(int) * 4)) 

printf ("Error writing file\n")/ 
return (1); 

) 

/* Now write out the kref vector */ 

bwritten » write (fhndl, (char *)KrefTab, sizeof (float) *Kldx) ; 
if (bwritten 1- (aiaeof (float) *KIdx) ) 

printf ("Error writing file\n M ); 
return (l); 

> 

/* output bin count */ 

printf ("Total Kref bin count :%ld\n",binont) ; 



Detection algorithm for above histogram 

*/ 

Function: 

* DoBox 

* De a exit ion: 

* Process the box. 

* Usage: 

* DoBox (x, y) 

* Inputs; 

* x - int : the x coordinate of the candidate pixel 

* y - int : the y coordinate of the candidate pixel 
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* Outputs: 

* Hone 
V 

static void DoBox (int x, int y) 
int tx, ty; 

double diffK, diffL, diffK; 
double kreflo,krefhi,krefavg; 
/* int tap; */ 
double mindiff; 

Pixel *midpxl m &8oanLine[y][x); 
Pixel *pxl; 

/* Averaae the values for this pixel */ 
AveragePixel (x, y) ; 

/* See if ve need to do this pixel */ 
if (midpxl->avghia > 2000.0) 
return; 

/* 

* Calculate the min difference value (this is calculated by using 

• twice the expected noise as the difference value). 
*/ 

mindiff - (10000.0/ (100. 0+midpxl->avghia) ) ; 

/* Nov loop through the pixels doing the box */ 
for (ty - y - BORDER; ty <- (y + BORDER) ; ty++) 

/• Get the pixel */ 

pxl - £ScanLine(ty][x • border]; 

/* Loop through the x */ 

for (tx • x - BORDER; tx (x + border) ; tx++, pxl++) 

/* See if ve need to look at this pixel (edge* are no-no c) */ 
If (pxl->sobel) 
continue; 

/* Average this sucker */ 
AveragePixel (tx, ty); 

/* Nov difference the Hi AIRS */ 
diffH « midpxl->avghia - pxl->avghia; 

/* Now threshold it */ 
if (diffH < mindiff) 
continue; 

/* Nov difference the Lo AIRS */ 
diffL = uidpxl->avgloa - pxl->avgloa; 
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/* Now threshold it */ 
ir ((diffL < aindiff) || (diffL — 0.0)) 
continue; 

kref lo-LookupKref (pxl->avghia , pxl->avgloa) ; 
krefhi-LookupKref (midpxl->avghia,midpxl->avgloa) ; 

dlffK - diffH/diffL; 

/* Key lookup algorithm 
* Histogram generation algorithm has been fit to this ratio 



kref avg- (kref lo* • 8 ) + ( . 7 *krefhi) ; 

/* See if ve need to histogram this point */ 
if (fdiffK < (krefavg+(KinThreshold))) 

f | (dif fK > (kref avg+ (MaxThreshold) ) ) ) 

continue; 

midpxl«->histval++ ; 

> 

) 

if (maxhit<midpxl->histval) 

ttaxhit»aidpxl->histval ; 
if (aidpxl->histval > fomThreeh) 

foa 4« (midpxl->hi8tval - f o&Thresh) ; 

if ( (midpxl->histval > 0) &6 (midpxl->hietval <200)) 
histpix[midpxl->histval] ++; 



> 
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